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In  this  report  ve  first  determine  representations  for  the  Anger- 
Weber  functions  «Jy(ax)  and  IE^(ax)  in  series  of  symmetric  Jacobi  poly¬ 
nomials.  (These  Include  LegSndre  and  Chebyshev  polynomials  as  special 
cases.)  If  v  is  an,?  integer,  these  become  expansions  for  the  Bessel 
function  of  tne  first  kind,  since  JFQ(ax)  »  Jn(ax)  .  Next,  corresponding 

representations  are  found  for  (ax)“*J^(ax)  .  Convenient  error  >ynpdf  are 
obtained  for  the  Chebyshev  cases  of  the  above  expansions . 

In  the  final  section  of  the  report  ve  determine  the  similar  type 
expansions  for  the  Bessel  functions  Yn(ax)  and  K^ax)  . 
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1.  Symmetric  Jacobi  Expansions  of  Anger-Y/eber  Functions 


A  function  f(x)  satisfying  certain  conditions  (for  these  consult 
[l]*)  may  he  expanded  in  the  series 


f(x)  »  ,  -l*x£l  ,  oC>-l  ,  (1.1) 

n«*0 


where  P^^^x)  is  called  the  symmetric  Jacobi  polynomial  of  degree  n  . 
For  our  present  purposes  we  shall  use  a  definition  given  in  [2]  : 


2nnJP^0C>(x)  -  (-Pd-x2)-**1  [(l-x2)0^11]  .  (1.2) 


Also 


f(x)(l-5i?fp(<^Gt')(x)dx  , 


(1.3) 


2goC  r(n+g:+l)  . 


n!  ^n+oC+  (n+cc+1)^ 


■^f '  *  ’ 1  •  (1-4) 


Using  the  representation  (1.2)  in  (1.3)  and  noticing  that  all 
derivatives  of  (l-x2)oC+n  up  to  and  including  the  (n-l)st  vanish  at 
x  «  tl  ,  we  integrate  (1.3)  n  times  by  parts  to  get: 


Cn  -  (2nnJhn)'1^  f(n)(x)(l-x2)a+ndx  .  (1.5) 


*  Numbers  in  brackets  refer  to  bibliography  at  end  of  report. 

This  report  was  released  by  the  author  on  13  December  1960  for  publication 
as  a  WADD  Technical  Report. 
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Consider  the  integral  definition  of  the  Anger-Web er  functions  [3] 

j^(ax)  +  il^ax)  -  Tf1  J  e*  8in  «  f(x)  ,  (1.6) 


When  y  is  an  ihteger  Jv(ax)  coincides  with  the  Bessel  function  of  the 
first  kind  Jv(ax)  [4]. 

Now  differentiate  (1.6)  n  times  under  the  Integral  sign,  sub¬ 
stitute  the  result  in  (1.5)  and  interchange  the  order  of  integration  (which 
is,  of  course,  permissible).  The  inner  integral  is  known  [5]  and  after 
evaluating  it  we  have 


Cn  ■  ( 


•i )“(*»! k  2  y1  J  e1  | 


fa  sin 

l  2 


(~£) 


X  J  1  (a  sin 
n+oc+  — 

2 


(1.7) 


Use  the  power  series  expansion  for  the  Bessel  function  in  (1.7)  and  inte¬ 
grate  tens  by  term  to  get 

cn  -  (-*)“  [cos  2Z  +  i  sin  2^  J  A^V,  oC,a)  ,  (1.6) 


where 


A_  ■ 


P(f +  +l)  r(|  -  f  +1) 


(1.9) 


and  Rn  is  conveniently  described  in  hypergeometric  notation  [6]  as 
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Rn(V,flC,a)  -  2F3f|  ♦  |  ,  |  +lja+n+  |  ,  |  +  g  +1,  §  -  |  +U-  ~1  (X.10) 

►  • 


Equating  real  and  laaglnary  parts  of  (1*6)  and  (1.2)  we  get 


JL(ax)  »  X  Vb^(x)  >  -X^x^X  , 
n»0 

(1.11) 

Ev(ax)  -  BbP|“»*)(x)  ,  -lSxSl  , 

n»0 

(1.12) 

where 

An  -  ac^aJAi^)  > 

(1.15) 

and 

Bn  ■* 

^.nRn(V,oC,a)yn(V)  ,  , 

n 

(-)^  com  2i2£  ,  n  even  , 

(1.1*) 

A,M  -  - 

n-1 

(-)  ®  sin  ,  n  odd  j 

■  n 

(-J2  sin  2L2L  ,  n  even  , 
o 

(1.15) 

vn<v>  ■  ' 

c 

n+1 

(-)  2  cos  2LZ  ,  n  odd  . 

1  2 

(1.16) 
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When  cc  ■  -  i  ' 


1 

r C1*  2  )  (ni1C* )  *n(*)  »  »  -  4#2. . .  ,  (1.17 ) 


vfaere  Tn(x)  is  the  Chebyebev  polyncsAal  of  the  first  kind  of  degree  n  . 
Also  for  this  value  of  oe ,  simplifies  to  the  products  of  two  Bessel 

functions  [7].  With  a  -  -  i  ,  then  (l.U)-(l,14)  become 


■  Z  <**»(*>  ,  -l£x£l  ,  (1.16) 

n*0 


l£v(*x)  -  £  DnTn(x),  -lSxSl  ,  (1.19) 

n-0 


where 


°n  “  e°J£lv(|)  Jn^(|)  ,  (1.20) 

2  2 


^  *  €nJn+v  (f  )  Jn-»(|)  Vn(v)  ,  (1.21) 

2  2 

and  6n  -  [l,n»0j2#n>0]  . 


Tor  Integral  y  we  have  the  expansions 

•w**>  •  Z  enVn(|)  Jk-n(|)  Wx>  .  -ISxSl  ,  (1.S2) 
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^2k+l(ax)  "  2!  2Jk+n+l  (  ”)  ^k-n  (  f)  ^2n+l(x)  *  (1*25') 


n*0 


and  k  -  0,1,2 
Since 


•  •  •  •  • 


VTfl 

Jv(i*)  -  e  2  yt) 


(1.24) 


where  Iv(z)  is  the  modified  Bessel  function  of  the  first  kind  [a],  we  may 
replace  a  by  la  in  (1.22)  and  (1.23)  to  get  expansions  for  Igk(ax) 
and  l2fc+l(ax)  . 

*  .  ,v, 

,  It  is  important  to  note  that,  although  the  above  expansions  are 
valid  only  for  x  real  and  |x|£l  ,  (1.6)  is  entire  in  a  and  V  ,  and 
hence  a  may  be  chosen  arbitrarily  to  yield  expansions  valid  everywhere  In 
the  finite  complex  plane. 

The  expansions  (l.ll),  (1.12),  (1.1B),  (1.19),  (1.22)  and  (1.23) 
are  quite  rapidly  convergent,  particularly  in  the  Chebyshev  cases,  'Sat 
consequently  the  last  four  axpsuislans  are  eminently  suitable  for  *tfe  4* 
digital  computers.  Such  series  are  usually  truncated  and  rearranged  IBS 
powers  of  x  ,  although  this  is  not  necessary  since  Tn(x)  satisfies  simple 
recursion  relationships  [9]*. 

Let  )eH|  denote  the  maxi  mini  error  incurred  for  -l£x£l  by 
taking  Just  N  terms  of  any  of  the  expansions  (1.22)  and  (1.23).  Using 
[is]  and  the  inequalities 


|Tn(x)|Sl,  -lSxSl  ,  (1.25) 


(m*-n)J  ^minJmax  [m0,^]  ,  m,n  integers  not  both  sero  ,  (1.26) 


*  The  Bessel  functions  required  to  compute  the  coefficients  in  our  expan* 
a  ions  can  be  systematically  generated  on  electronic  computers  with  the 
aid  of  techniques  discussed  in  [lO,  11,  12].  There  are  numerous  tables 
available  for  hand  calculations. 
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v*  nay  derive  error  bound*  for  M>k  .  For  (1.22), 

|a|2N#xpM  (§)  I 


l€*l* 


24»-i(»fk)j  (H-k):  lo  C  |  TJ&jS  | )  (1#87  * 


end  for  (l.23) 


|c  ,>l"“*M(g)l  x  |  . 

*  24*fi(j|fk+i)j(B.k)!  0  I  2  4  (IWc+l)(*-k) 


JVar  the  cases  I gk(ax)  and  Ig^+iCax),  respectively,  exp  M(|)l 
(1.27)  end  (l.ae)  is  replaced  by  exp|2,fi  | 

2.  Expansions  of  Bessel  Functions  of  the  First  yina 
of  Ronlntegral  order 


Results  in  the  previous  section  gave  syneetrlc  Jacobi  poXjfo&iK&fe; 
expansions  for  J^(ax)  and  I^(ax)  for  Integral  y  ,  When  y  is  »m- 
integral,  these  functions  are  no  longer  entire  functions  of  x  ,  it  is 
convenient  to  study  the  entire  function 

r(v+l)(ax/2rV  Jv(ax)  -  qFj^^+1;- .  (2.1) 


we  will  derive  expansions  for  the  above  gF^  .  Corrsepoading  expansion* 
for  r(V+l)(ax/2)“V  Iy(ax)  then  follow,  as  before,  fTo«  (1.24). 

We  shall  first  need  the  following  identity. 

i  ■  i 

-  ( %")  cry1*?"  * }  p 
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where  the  generalised  Jacobi  polynomial  of  degree  n  ,  pjf^^x)  amy  be 
defined  as 

P^<*)  -  ("£*)  2h  n+Of +1; «+l;  |  -  |  x]  .  (8.3) 

The  following  quadratic  transformation  of  the  Gaussian  hypergeo- 
metric  function  is  given  in  [1*] : 

oFi  2a,2b;a+b+  i  j  !  -  i 

*  L  2  2  2 
1 

Let  (1-z)2  ■  x  ,  a  «  -n/2  ,  b  ■  n/2  +ce+  i  ,  and  c  »  oC+l*a  +  b  +  i. 
in  (2.4).  Then 

2F1  |^"njn+2<X-*-lj  a+1;  i  -  i  xj  *  2yl  ["  §  '  2  +CC+  2  '  Gf+lil-x2  J  ,(8.5) 

or,  if  we  refer  to  (2.3)  and  replace  n  by  2n  (2.5)  beccoms 

(x )  -  )  2F3.  [’n>  **  a+  |  >  «+lil-x2  ]  (2.6) 

and  this  is  simply  (2.2). 

If 

o*i jWj  *(^)1  -  f  2^(«)  *  (g-7) 

L  w  J  n-0 

then 

An-hJ1J  oFi[y.l;>(^)](l-sf(l.s)  2P^  '  2/(ej  #  (2.8) 


1 

(l-*)2 


J  “  2F1  j*  i*J  •  (2.4) 
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or,  equivalently, 


Aft  -  (h^i  r1  ^  <?1  [v>lj  -^-*]  (l-*)a+n(l+e)  *  de  ,  (2.9) 


(h^T~(g)a/1  Af^li^(l**)](^d)a+,l(l+»)  2  d*.  ,  (2.10) 


n-  i 


and 


*+ 1 

2  2  (ml). 


(*■*«*  i)  (m  £) 


(2.11) 


Nov  expand  the  hypergecroetrlc  function  in  (2.10)  in  its  power  series  and 
integrate  tens  by  term  to  get 


*nm 


- '  - .  jPo  [nt  i  J  »+m-l,2n +  c£+  £  ;  xl  •  (2.12) 

(m«.  i)^!^  1-2  2  J 


Let  z  ■  (2x?-l;  in  (2.7)  and  X  «*  -a2/4  >  and  use  (2.2).  Then 


jji*)  .(jif  § 


(2.1S) 


n-0 


where 


*n 


iifiar'2n 


*/7r(2n+2£Z+l), 


»£! - 1Fpfn+  i  ;  V+n+l,2n+rf+  5;  -  sf|.(2. 

L  ('mi')  1 2L  2  2  ‘  J 

2l1^  2'>)+i 
2 


M) 
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For  the  Chebyshev  case  c c 


1  and 

2 


M“>  .(H)V  £  Vsa(x)  , 

n-o 


-Kxll  , 


where 


ni  f(J>+n+l)  1  ® 


1  j  V+n+l,2n+l; 

2 


(2*15) 


(2.16) 


Notice  that  when  V  «  -  i  ,  (2.14)  simplifies.  Also,  since 

2 

_  1 

J  1  (ax)  2  cos  (ax)  (2.17) 

*2 


we  infer  the  expansion  - 


cos  (ax)  «  £  cnpln,tr^x) 


(2.18) 


n*-0 


where 


■Cr 


i  -«+  i 

(-)n7f22  (^2n+a:+  i)(2n+cr+l)a 

~T~ 

a  ^ 


2n+  et+  - 
2 


(a)  ,  (2.19) 


a  formula  which  can  be  derived  in  a  number  of  different  ways. 

Another  Interesting  and  quite  rapidly  convergent  expression  for 
the  coefficients  in  (2.13)  can  be  derived  by  substituting  the  expansion  (see 
Watson  [l3,  p.  142]) 
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2  \k 


00  (§)  W*)*1-*  >' 


k=0 


kJ' 


(2.20) 


for  f(x)  in  (1.3).  Employing  a  known  integral  formula  {l5j,  we  have 


2V“2Cr’l'vnt(4n+2a:+1)(2n+ar+l)ar  «  $n,k(§)  JD+k^a) 

®n  = - - - 2_ 


b?  (2n+l) 


M)  (cc+k+l) 


(2.21) 


where 


$n,k  "  3F2(-2n,2n+2<r+l,  ar+k+1;  oC+l,2k+2ar+2;l)  .  (2.22) 


If  the  expansion  (2.13)  is  truncated  after  N  terms,  we  can 
obtain  rough  error  estimates  for  large  N  by  assuming  that  the  expansion 
converges  rapidly  enough  so  that  the  first  neglected  term  is  approximately 
the  error  incurred.  Using  the  relation 


(2.23) 


we  see  that  the  jF2  in  (2.14)  may  be  replaced  by  an  exponential  tern. 

Also 

lp!?a)(*)l  *  CiT)’  "1*x*1  '  <S*24) 


so  for  (2.13)  we  have  roughly 


a2 

‘ 

e"  8N 

P(2w-2  or+i) 

24Nf2hr+*>  p^2W+oP+  i)r(R*-V+l)rV+l) 


-l^x^l 


(2.25) 
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3*  Expansions  of  Bessel  Functions  of  the  Second  Kind 


The  Bessel  function  and  modified  Bessel  function  of  the  second 
kind  are  denoted  by  Yv(z)  and  Kjz)  ,  respectively,  and  a  treatment  of 
them  can  be  found  in  [2,  v.  2,  Ch.  VII  ].*  If  p  ls  nonintegral,  then 

-  [sin(V7T)]  1  |jy(z)  cos(vr)  -  J„v(z)|  ,  (3.1) 


Mz)  -  ^[ainOJr)]  _1  jl_„(z)  -  IJt)J 


M) 


so  for  such  values  of  p  expansions  for  the  functions  follow  directly 
the  results  of  section  2. 

If  V  is  an  integer,  it  can  be  shown  that 


V“>  =|r  J*(«*»og(jK)  +  Vl<“>  *  V°k(ax)  ' 


Kk(ax)  »  (-)k*\(m)lce(Si).  X  i^^dax)  ♦  G^lax)  ,  (5.4) 


where 


,Nk.i(ax)  = 


_  K-J.  , 

i-  V  z  Cf  )■ 

m=n 


k>o 

m  * 


k  «  0 


(3.5) 


*  This  reference  calls  K  (z)  the  modified  Bessel  function  of  the  thitti 
kind.  ■  1  - 
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I 


and 


Gjax)  -  V  (-)nf^k+2,B  Lv<(k-fnH-lK  ^(**1)3 

'  al  (k+m)! 


(3.6) 


We  assume  the  value  of  log  ^ Is  known.  Then,  since  expan¬ 
sions  for  J^(ax)  and  I^(ax)  were  found  in  section  1,  and  since 
®k-l(ax)  aiaply  a  polynomial  in  l/(ax)  ,  we  need  expand  only  the  entire 
part  of  (3.3),  i.e.,  Gjj(ax)  ,  in  synmetric  Jacobi  polynomials. 

Using  the  representation  (3.6)  as  f(x)  in  formula  (1.5),  a 
straightforward  derivation  gives  the  series 


Gk(ax)  «  £  AnPfo^x)  ,  -1<x<1 

n=0 


(3.7) 


where 


An  = 


[(-^(-^(n+of+l)^  (n+<*:+ |) 


(-^(-k-Pm^ 


3n+2«f+l 


m-0 


n+tt+1 


v  /  a  \k+2m  [y(k+nn-l)+  vKnrfl)] 
*  ^2'  m:(k+m)! 


(3.8) 


4.  Conclusion 


The  symmetric  Jacobi  expansions  given  in  this  paper  converge, 


in  general,  most  rapidly  when  oc 


i.e.,  for  the  Chebyshev  case. 
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Another  advantage  gained  when  oc  =  -  i  is  that  the  hypergeometric  func- 

2 

tions  occurring  in  the  coefficients  of  the  expansions  in  several  cases 
simplify  considerably. 

5y  using  the  same  techniques  expounded  in  this  paper,  expansions 
of  Bessel  functions  in  Iaguerre,  and  hence  by  a  simple  transformation, 
Hermite  polynomials  can  be  obtained,  but  the  convergence  of  these  expansions 
is  vastly  inferior  to  that  of  those  derived  here. 

Since  our  expansions  converge  much  more  rapidly  than  the  corres¬ 
ponding  Taylor's  series  representations,  they  will  be  of  great  utility  as 
subroutines  for  evaluating  Bessel  functions  on  digital  computers.  We 
emphasize  the  need  far  calculating  in  the  near  future  the  coefficients  of 
these  expansions,  particularly  for  the  Chebyshev  cases. 
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